The most probable wave function of a single free moving particle 
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We develop the most probable wave functions for a single free quantum particle given its mo- 
mentum and energy by imposing its quantum probability density to maximize Shannon information 
entropy. We show that there is a class of solutions in which the quantum probability density is 
self-trapped with finite-size spatial support, uniformly moving hence keeping its form unchanged. 
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Introduction: Madelung Fluid — Since its final formu- 
lation in term of Schrodinger wave mechanics, quantum 
mechanics has claimed to have never failed any conceiv- 
able experimental test Yet, despite this overwhelming 
fact, quantum mechanics could not hide its embarrassing 
face that it can not give an unambiguous answer to a 
very simple yet fundamental question: what is the wave 
function of a single free particle moving linearly (trans- 
lationally) with finite velocity. To describe a single free 
particle, one can not simply use the stationary plane wave 
for the reason that it extends to the whole space so that 
formally unnormalizable, besides is physically in contrast 
to our observation that a particle is spatially very well 
localized. In his correspondence to Einstein, Born sug- 
gested to use coherent wave packet [2] , motivated by the 
fact that a coherent state is localized both in momen- 
tum and position space. Yet, coherent wave packet suf- 
fers from the lack of being not stationary to be evolved 
by Schrodinger equation. This is again is in contrast to 
our intuition about a particle which should have a stable 
structure, if left unperturbed. 

On the other hand, in our previous work [3], working 
with two dimensional Madelung fluid dynamics Q whose 
irrotational motion reduces into the Schrodinger equation 
for a single free particle, we have developed a class of the 
most probable wave functions for a single free particle 
with a given quantum energy. This is done by imposing 
its quantum probability density to maximize Shannon 
information entropy. We showed that there is a class of 
solutions in which the wave function is self-trapped with 
finite-size spatial support, spinning around its center yet 
stationary. In this present paper, we shall develop yet 
another class of stationary solutions in which the wave 
function is still self-trapped with finite-size support and 
moving linearly with uniform (thus definite) velocity. We 
then end to give an interpretation to the wave function 
as describing a single free particle moving with finite and 
definite velocity. 

Let us consider the Madelung fluid dynamics of a sin- 
gle free particle with mass m. Since in the present paper 
we are only interested in the translational motion of the 
particle, it is then sufficient to consider a spatially one 
dimensional case. First, the state of Madelung fluid in 
configuration space q is determined by a pair of fields, 
{p(q),v(q)}. p(q) is called as quantum probability den- 
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FIG. 1: The profile of a self-trapped quantum probability den- 
sity p(q) (solid line) and the corresponding quantum potential 
U(q) (dashed line) that it generates. The quantum potential 
is shifted down so that its global minimum is equal to zero. 
See text for detail. 



sity thus assumed to be always normalized. Whereas v(q) 
is a velocity vector field. Their evolution in time t is then 
governed by the following coupled differential equations: 



m ~oTt = ~ dqU ' dtp ' 



d q (vp)=0. 



(1) 



Here, U (q) is called as quantum potential determined by 
the quantum probability density as 



U(q) 



h 2 d 2 q R 
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(2) 



where R(q) 



A/2 



(q) is quantum amplitude. Notice that 



for one dimensional system, one can always derive the 
velocity field from a scalar field S(q) as v(q) = d q S/m. 
Using this to define a complex-valued wave function as 
ip(q) = R(q) exp(iS(q)/h), Madelung fluid dynamics of 
equation (pQ) is then equal to the Schrodinger equation 
for a single free particle: 



ihd t ^(q;t) 
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(3) 



In this regard, S(q) is called as quantum phase. 
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The most probables elf -trapped quantum probability 
density — Next, let us consider a class of quantum prob- 
ability densities which maximizes Shannon information 
entropy: H[p] = — J dq p(q) lnp(g) [5[, given its average 
quantum potential U = J dq U(q)p(q). This is the so- 
called maximum entropy principle [6] . It has been argued 
as the only way to infer from an incomplete information 
which does not lead to logical inconsistency pj. Hence, it 
will give us the most probable quantum probability den- 
sity with average quantum potential U. This inference 
problem can be directly solved to give [8(]: 

p(q) = ^exp(-U(q)/T), (4) 

where T is the Lagrange constant below chosen to be non- 
negative, and Z(T) is a normalization factor. Combining 
with equation ([2j), equation (|4]) comprises a differential 
equation for p(q) or U (<?), subjected to the condition that 
p(q) must be normalized. In term of U(q), one has Q 

d 2 q u = ±{d q u? + ^fu. (5) 

Figure [U shows the solution of equation ([5]) with the 
boundary conditions: 17(0) = 1, d q U(0) = 0, for T = 1. 
All numerical results in this paper are obtained by set- 
ting h = m = 1. One can see that the quantum proba- 
bility density is being trapped by its own self-generated 
quantum potential. The self-trapping property is valid 
for any positive value of £7(0) provided that d q U(0) = 0. 
This can be understood as follows. First, at the origin, 
q = 0, we have d%U(0) = (4raT / h 2 )U (0) > 0, hence U(q) 
is convex at q = 0. Since d q U(0) = 0, then at points 
nearby q = 0, one concludes that U(q) is positive defi- 
nite. Furthermore, since the first term on the right hand 
side of equation ([5]) is always non- negative, then at these 
points one has d 2 U > 0. This reasoning can be extended 
to the whole space so that U(q) is everywhere convex 
and positive definite. Using this fact in equation (jH), one 
concludes that p(q) is trapped by its own U(q). 

Next, from the definite positivity of U(q), then the 
definition of quantum potential of equation ([2|) gives us 

d 2 q R < 0. (6) 

R(q) is thus concave everywhere. Since R(q) is finite 
and possesses symmetry R(—q) = R(q), then R(q) must 
cross the q— axis at finite values of q = ±L m where 
U(±L m ) = oo. See Fig. [TJ Hence, the self-trapped quan- 
tum probability density developed in previous paragraph 
possesses only a finite range of spatial support given by 
the interval: M = [— L m , L m ]. 

In Fig. [2j we plot (a) L m , and (b) the average of 
the quantum potential U = J dqll p, of our self-trapped 
quantum probability density against the variation of the 
parameter T, by solving equation ([5]) with fixed boundary 
conditions: d q U(0) = and £7(0) = 1. Two asymptotic 
limits are of great interest. First, for the case when T 
takes an infinite value, T — > oo, L m is approaching zero, 




FIG. 2: (a) Half length L m of the support of self-trapped 
quantum probability density, and (b) its average quantum 
potential U plotted against T. 

liniT^oo L m = 0. One thus concludes that the quantum 
probability density is approaching a delta function for 
sufficiently large value of T. Further, one can also see 
that in the limit T — > oo, the average quantum potential 
is infinite, liniT^oo U — > oo. Based on this last fact, we 
shall show later that the limit T — > oo is not relevant 
physically. 

On the contrary, for vanishing value of T, we see from 
Fig. [2] that both L m and U are approaching some fi- 
nite values, hniT^o = To and liniT^o U = Uq. This 
hints us that the quantum probability density and thus 
its corresponding quantum potential are also converging 
toward some functions 

limU(q;T) = U (q), limp(q;T)=p (q). (7) 

Let us discuss this asymptotic situation in more detail. 
To do this, in Fig. [3] we plot the profile of quantum prob- 
ability density and its corresponding quantum potential 
for several small values of T, T = 1, 0.5, 0.05, by solving 
equation ([5]) with fixed boundary conditions: d q U(0) = 
and £7(0) = 1. Let us first pay our attention to the pro- 
file of quantum potential U(q;T). One can see that as 
T decreases, the quantum potential inside the support is 
getting flatterer before becoming infinite at the bound- 
ary points, q = ±L m (T). One can thus guess that at 
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FIG. 3: The quantum probability density (solid line) and 
its corresponding quantum potential (dashed line) for several 
small values of T obtained by solving equation J5}. We also 
plot the analytical solution for po(q) at T = 0, assuming that 
Uo(q) takes the form of a box with infinite wall. See text for 
detail. 

vanishing T, T = 0, the quantum potential is converging 
toward a box of length 2Lq with perfectly flat bottom 
and an infinite wall at the boundary points, q = ±Lq. 
See Fig. El 

Based on the above guess, let us proceed to calculate 
the profile of the quantum probability density for van- 
ishing T, po(q). To do this, let us denote the assumed 
constant value of quantum potential as Uo(q) = U c , for 
— Lo < q < Lq. Since Uo(q) is flat inside the support, 
then one has Uo = dqpo(q)Uo(q) = U c . Hence, from 
the definition of quantum potential given in equation ([2]) , 
inside the support of the quantum probability density, 
one has 

ti 2 

— - — d R{) = U c Rq = UoRq, (8) 
zrn y 

1/2 

where Ro(q) = p Q (q) is the quantum amplitude at T = 
0. The above differential equation must be subjected to 
the boundary condition: Rq(±Lq) = 0. 
Next, solving equation JS|) one obtains 

Ro(q) = A cos(k q), (9) 

where A$ is a normalization constant and ko is related to 
the average quantum potential as 

k = ^2mU /h 2 . (10) 

The boundary condition imposes koLo = tt/2. In Fig. 02 
we plot the above obtained quantum probability density, 
po(q). One can see that as T is decreasing toward zero, 
p(q;T) obtained by solving the differential equation of 
(j5j) is indeed converging toward po(q) given in equation 
([9]). This confirms our initial guess that at T = the 
quantum potential takes a form of box with infinite wall 
at q = ±Lq. 



Stationary-moving wave packet — Let us now choose 
the pair of fields {po(q),vo(q)} as the initial state of the 
Madelung fluid. Here, po(q) = Rq(q) is given in equation 
(|9|) and vo(q) = v c is a uniform velocity field with non- 
vanishing value only inside the spatial support M. Then, 
since at t = the quantum potential is flat inside the sup- 
port, the quantum force is initially vanishing, d q U = 0. 
Hence, from the dynamical equation (pQ), the velocity field 
at infinitesimal lapse of time t = At will stay unchanged 
and keeps uniform. This in turn will shift the quantum 
probability density in space by Aq = v c At while keeps its 
spatial profile unchanged. Correspondingly, it will shift 
the support of the quantum probability density by the 
same amount: M. — > M/±t = [v c At — Lo,v c At + Lq\. 
The same thing occurs for the next infinitesimal time 
lapse and so on and so forth. Hence, at finite time lapse, 
the velocity field is kept uniform with constant value. 
This fact finally gives us 

dtp = dtp + v d q p = dtp + d q (v p) = 0, (11) 

where in the second equality we have used the continuity 
equation of (pQ) and the fact that the velocity field is uni- 
form so that d q vo = 0. One thus concludes that the quan- 
tum probability density is moving with a uniform velocity 
field v c , keeping its initial form remained unchanged. 

Hence, at time t, the quantum probability density is 
given by 

p(q; t) p (q - v c t) Aq cos 2 (k q - u t), (12) 

where c^o = kov c and q G [v c t — Lo,v c t + Lq\. One can 
check easily that equation ([T2]) indeed satisfies the con- 
tinuity equation of (pp). Moreover, the quantum phase at 
time t can be obtained by integrating v(q) = (l/m)d q S = 
v c to give 

Sfat)=mv c q + Z(t), (13) 

where is a function which depends at most only on 
t. We shall give the explicit form of £(t) later thereby 
revealing its physical meaning. The wave function we 
are searching thus takes the following form: 

V>(<?; t) = A cos (k (q-v c t)) exp Q (mv c q+£(t))^) , (14) 

where q G [v c t — L , v c + L ]. 

Next, let us calculate the quantum mechanical en- 
ergy of our self-trapped, uniformly moving-stationary 
wave function. Since quantum mechanical energy is 
conserved then it is sufficient to use the wave function 
at t = 0. Putting the wave function in polar form, 
ip(q) — Ro(q) exp(iS(q)/h), one has 

(E) = J L ° dq r(q)(-^d^(q) 

- l ^R, dq R 0dq S-^-RldlS). (15) 
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The first term on the right hand side is equal to 
the average quantum potential, Uq = J dqUopo = 
tfk 2 ,/ (2m). Next, defining kinetic energy density as 
Ko(q) = (ra/2)i>o(<z), then the second term is equal to 
the kinetic energy of the Madelung fluid: 

K = J dqK (q)p (q) = (16) 

Further, for a uniform velocity field, the last term is van- 
ishing, (l/m)d 2 S = d q vo = 0. Again for a uniform ve- 
locity field, since Ro(q) is an even function and d q Ro(q) 
is an odd function then the third term is also vanishing. 
Hence, in total, the quantum mechanical energy of the 
self-trapped wave function moving with a uniform veloc- 
ity field can be decomposed as 

W=«. + ft = ^ + =£ (17) 

Moreover, using similar argument as above, the average 
quantum momentum can also be calculated to give 

(p) = [ ° dq r(q)(- ihd q )^j(q) = mv c . (18) 

The average kinetic energy of the Madelung fluid and 
average quantum momentum are thus related as Kq = 
(p) 2 /(2m). This observation leads us to conclude that 
Uq must essentially be interpreted as an internal energy 
of the single particle. It is this energy that is missed to 
be taken into account if one uses the simple plane wave 
to represent a moving particle. 

Finally, putting the wave function of the form ([T4|) into 
the Schrodinger equation of (J3]), keeping in mind uoq = 
kovo and equation (jTTJ), one easily obtains dt£ = —(E), 
which can be integrated to give 

£(t) = -(E)t, (19) 



modulo to some constant. Putting everything into equa- 
tion (fT3]h one finally obtains 

1>(q; t) = A cos (k (q - ^i)) exp (i ({p)q - (E)t)), 

(20) 

where q G [(p)t/m — Lq, (p)t/m + Lq]. Hence, since ko 
is determined by fixing the values of (E) and (p), the 
stationary-moving wave function we developed in this 
paper can be interpreted as describing the most proba- 
ble wave functions for a single free moving particle with 
momentum (p) = mv c and total energy (E). 

Conclusion and Discussion — To conclude, we have de- 
veloped a class of the most probable wave functions for a 
single quantum particle moving with finite velocity. This 
is done by imposing the wave function to maximize Shan- 
non information entropy. We showed that there is a class 
of solutions in which the wave function is self-trapped by 
its own quantum potential such that it is localized with 
finite-size spatial support; translationally moving with 
uniform velocity thus keeping its form unchanged. We 
also showed that there is a new type of energy, referred 
to as internal energy, which is missed if one uses a plane 
wave to represent the free moving particle. 

The fact that the velocity field is uniform in space as- 
serts that the momentum of the wave packet is definite. 
Surprisingly, even in this case, the wave packet possesses 
only a finite uncertainty in position. This observation 
thus raises an interesting question concerning the physi- 
cal status of Heisenberg uncertainty principle. 
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